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ABSTRACT 

This  paper  provides  some  background  and  historical  review  of  how  microstructure-sensitive  finite 
element  simulations  can  play  a  role  in  understanding  effects  of  stress  amplitude,  R-ratio,  and 
microstructure  on  fatigue  crack  formation  and  early  growth  at  notches,  including  pores  and  non-metallic 
inclusions  for  Ti  alloys  and  Ni-base  superalloys.  Fatigue  indicator  parameters  (FIPs)  are  computed  over 
finite  volumes  that  relate  to  processes  of  fatigue  crack  formation  and  early  growth  at  the  scale  of 
individual  grains.  It  is  argued  that  both  coarse  scale  (uncracked,  mesoscale)  and  fine  scale  FIPs 
(computed  in  the  vicinity  of  cracks  in  single  grains  or  crystals)  can  be  related  to  the  cyclic  crack  tip 
displacement  that  serves  as  a  driving  force  for  crystallographic  fatigue  crack  growth,  and  that  the  fine 
scale  FIPs  correlate  directly  with  cyclic  crack  tip  displacement.  Scatter  in  HCF  and  VHCF  is 
computationally  assessed  using  multiple  statistical  volume  elements  and  the  distribution  of  FIPs  of 
extreme  value  character.  The  concepts  of  marked  correlation  functions  and  weighted  probability  density 
functions  are  reviewed  as  a  means  to  quantify  the  role  of  multiple  microstructure  attributes  that  couple  to 
enhance  the  extreme  value  FIPs  in  the  HCF  regime.  An  algorithm  for  estimation  of  the  cumulative 
probability  distribution  of  cycles  for  crack  formation  and  growth  from  notches  in  HCF  and  VHCF  is  also 
described. 

1.  Background  and  Objectives  of  Microstructure-Sensitive  Fatigue  Modeling 

Quantifying  variability  of  fatigue  life  in  the  HCF  and  VHCF  regimes  is  a  common  goal  to  support 
design  based  on  low  probability  of  failure  and  typically  requires  extensive  and  time  consuming 
experimental  studies.  We  primarily  consider  here  the  effects  of  microstructure  on  fatigue  variability  of  Ti 
alloys  and  Ni-base  superalloys,  advanced  alloys  used  in  aircraft  gas  turbine  engine  components.  The  goal 
of  microstructure-sensitive  polycrystal  plasticity  simulations  is  to  assist  in  reducing  the  number  of 
experiments  required  by  developing  computational  schemes  that  can  be  used  to  help  design,  analyze  and 
qualify  new  microstructures  that  are  more  fatigue  resistant  for  a  given  application.  Specifically,  the  intent 
is  to  characterize  the  influence  of  microstructure  on  fatigue  life  variability  under  constant  amplitude 
loading  conditions,  with  particular  attention  to  the  shape  of  the  low  cumulative  probability  tail  regions. 
An  added  benefit  is  a  direct  pathway  to  identifying  mesoscale  microstructure  attributes,  either  by 
themselves  or  via  coupling  amongst  multiple  attributes,  that  contribute  to  “hot  spot”  neighborhoods  with 
high  probability  form  fatigue  cracks. 

The  focus  of  the  present  state-of-the-art  is  on  the  role  played  by  mesoscale  arrangement  of  grains  and 
phases,  encompassing  various  spatial  statistics  including  orientation  and  disorientation  distributions  [1]. 
While  variability  also  may  arise  from  the  nanoscopic  scale  in  association  with  crack  nucleation 
mechanisms,  it  is  useful  to  consider  that  much  of  the  observed  variability  in  polycrystalline  systems  arises 
from  microstructure  at  scales  well  above  1  |Lim,  considered  here  as  scale  representative  of  embryonic  stage 
of  crack  formation.  The  origins  of  these  kinds  of  mesoscopic  simulations  of  heterogeneous  materials  to 
understand  the  role  of  microstructure  and  related  stochastic/probabilistic  aspects  with  applications  to 
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various  responses  and/or  properties  can  be  traced  back  to  the  early  to  mid-1980s  with  the  advent  of 
computational  micromechanics.  For  example,  Asaro  [2],  Asaro  and  Needleman  [3],  McHugh  et  al.  [4], 
and  Deve  et  al.  [5]  developed  early  polycrystalline  finite  element  models  for  shear  localization  in 
polycrystalline  microstructures  with  ductility  and  strength  in  mind.  Needleman  and  colleagues  [6-7] 
considered  the  growth  of  grain  boundary  cavities  and  interactions  of  fibers  with  matrix  material  using 
finite  element  methods  in  the  process  of  nucleation  and  growth  of  voids.  Needleman  [8]  extended  these 
approaches  to  address  interfacial  failure  in  finite  element  models  of  idealized  or  actual  microstructure 
representations  using  cohesive  zone  methods.  By  the  early  1990s  considerable  efforts  were  underway  in 
computational  mesoscopic  modeling  of  microstructure-property  relations  directed  towards  various  classes 
of  heterogeneous  materials.  Extending  earlier  works  of  Asaro,  Needleman  and  colleagues,  Zikry  [9-11] 
considered  the  effects  of  crystallite  orientation  distribution  and  grain  boundary  network  on  failure  modes 
in  large  scale  computational  crystal  plasticity  using  finite  element  methods. 

The  integration  of  computational  micromechanics  with  probabilistic  failure  descriptions  was  a  natural 
outcome  of  the  rapid  development  of  the  field  of  computational  micromechanics  in  the  1980s.  Haddad 
(1990)  considered  the  effect  of  microstructure  on  distribution  of  slip  in  poly  crystals.  The  work  of  Chamis 
and  colleagues  [13]  at  NASA  integrated  the  notions  of  computational  micromechanics  simulation  of 
heterogeneous  materials  (composites)  with  probabilistic  analyses  to  advocate  low  probability  of  failure 
design  methods  for  different  failure  modes.  Duva  et  al.  [14]  and  Ruggieri  and  Dodds  [15]  employed 
computational  simulations  of  heterogeneous  materials  and  cracked  to  support  probabilistic  computational 
micromechanics  characterization  of  damage  evolution. 

Appealing  to  much  earlier  work  on  stochastic  deformation  and  failure  of  random  heterogeneous 
materials  (cf.  Krajcinovic  [16]),  and  in  recognition  of  the  expense  of  finite  element  methods,  lattice¬ 
spring  methods  were  used  by  Ostoja-Starzewski  et  al.  [17]  to  characterize  probability  distributions  of 
damage  evolution  in  random  heterogeneous  materials.  Zhou  and  Zhai  [18-19]  built  on  the  earlier  works 
of  Needleman  to  model  dynamic  fracture  of  realistic  microstructures  by  placing  cohesive  zone  elements 
between  all  elements  in  explicit  finite  element  simulations,  including  effects  of  distribution  of  interface 
cohesive  strength  and  energy.  Variability  of  response  (e.g.,  crack  growth  history)  is  a  natural  outcome  of 
such  simulations. 

Early  works  in  fatigue  of  heterogeneous  materials  were  typically  focused  on  composite 
micromechanics  [20-23].  Even  earlier  works  in  mesoscopic  simulations  of  metallic  polycrystals  for 
purposes  of  characterizing  distributions  of  driving  forces  for  fatigue  crack  formation  can  be  traced  to 
early  works  of  Pro  van  [24-25],  Cox  and  Morris  [26]  and  Sakai  et  al.  [27].  The  interaction  of  small  cracks 
with  microstructure  and  related  mesoscopic  models  for  crack  growth  in  random  microstructures  had  been 
understood  and  modeled  in  the  1980s  (cf.  Monte  Carlo  simulations  by  Tanaka  and  Akiniwa  [28]  with 
random  microstructure  that  correctly  predict  the  magnitude  of  scatter  of  the  small  crack  da/dN  using  a 
simple  da/dN  =  A(ACTOD)m  relation,  where  ACTOD  is  the  range  of  crack  tip  opening  displacement). 
McDowell  [29]  reviewed  the  status  of  microstructure  sensitivity  to  formation  and  growth  of  small  cracks 
in  HCF.  Gall  et  al.  [30]  used  finite  element  simulations  to  study  the  effects  of  crystallographic  orientation 
on  variability  of  microstructurally  small  crack  growth  rate  in  crystals  using  CTOD  concepts. 

The  mesoscopic  modeling  of  distributed  crack  formation  within  grains  of  a  polycrystalline  ensemble 
by  Hoshide  and  Socie  [31]  that  employed  slip  band  cracking  models  and  crack  growth  relations  was  a 
precedent  in  motivating  the  later  work  in  the  field;  all  essential  elements  of  crack  nucleation  in  slip  bands 
in  grains  and  early  growth  were  treated  for  metallic  polycrystals  with  the  approximation  of  Sachs 
(isostress)  among  grains.  McDowell  [32]  discussed  the  importance  of  such  mesoscopic  simulations  to 
understand  fatigue  phenomena.  This  work  clarified  that  for  fatigue  the  notion  of  a  statistically 
homogeneous  representative  volume  element  (RVE)  is  problematic;  many  practical  cases  it  is  too  large  to 
be  relevant  to  actual  specimens  or  structures,  particularly  in  the  HCF  regime.  This  had  been  stressed 
earlier  in  the  work  of  Lacy  et  al.  [33]  which  showed  that  the  RVE  size  for  elastic  stiffness  of  a  field  of 
microcracks  in  a  brittle  material  (e.g.,  polycrystalline  ceramics)  is  much  too  small  to  be  considered  as  a 
RVE  for  the  evolution  of  damage.  According,  it  is  necessary  to  build  up  statistics  based  on  multiple 
computational  realizations  [34]  of  smaller  statistical  volume  elements  (SVEs)  to  arrive  at  a  proper 
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statistical  distribution  of  responses,  or  to  employ  approaches  that  weight  input  from  SVEs  to  achieve 
much  more  efficient  characterization  of  ensemble  statistics  [35-36].  In  other  words,  it  is  not  appropriate 
to  select  a  size  for  a  computational  volume  of  a  heterogeneous  material  by  considering  convergence  of  the 
average  stress-strain  behavior  or  elastic  stiffness  if  fatigue  crack  formation  and  growth  are  the  relevant 
responses.  The  notion  of  a  RVE  implies  that  even  the  higher  order  moments  of  the  statistical  distribution 
of  fatigue  driving  forces  are  fully  captured,  not  just  the  lowest  order,  the  point  made  earlier  in  reference  to 
evolving  damage  [33]  and  fatigue  [34]. 

The  PhD  thesis  of  Bennett  [37]  and  derived  papers  by  Bennett  and  McDowell  [38-41]  used  finite 
element  models  with  crystal  plasticity  to  compute  distributions  of  slip  and  projected  distributions  of 
fatigue  cracks  in  polycrystals  based  on  multiple  realizations  of  a  nominal  microstructure.  In  addition, 
they  computed  the  cyclic  crack  tip  displacement  range  (both  opening  and  sliding  components)  as  a 
function  of  crack  length  relative  to  grain  size  from  surfaces  in  polycrystals.  In  the  DoE  USCAR  program 
from  1995-2000,  McDowell  and  colleagues  [42-46]  employed  finite  element  modeling  at  the  scale  of 
microstructure  to  model  formation  and  early  growth  of  interdendritic  cracks  in  A356-T6  casting  alloys, 
culminating  in  the  formulation  of  a  probabilistic  framework  for  fatigue  response  based  on  various 
distributions  of  microstructure  attributes  at  five  different  length  scales  of  attributes.  During  the  same  time 
period,  Tryon  and  Cruse  [47-49]  related  the  notion  of  mesoscopic  modeling  of  microstructure  to 
reliability  within  the  context  of  finite  element  simulations  of  polycrystalline  microstructures. 

The  past  decade  has  witnessed  a  significant  increase  in  attention  devoted  to  mesoscopic  modeling  of 
cyclic  deformation  processes  at  the  microstructure  scale.  In  addition  to  the  scheme  outlined  in  the  next 
section,  recent  works  in  this  regard  by  Dunne  and  colleagues  [50-55]  and  Ghosh  and  co-workers  [56-59] 
are  noted. 

2.  FIP  Methodology 

Microstructure-sensitive  crystal  plasticity  modeling  methods  have  rapidly  developed  within  the  past 
decade  [1,46,60-62].  The  basis  for  these  methods  is  consideration  of  effects  of  microstructure 
morphology  (size,  shape,  orientation/disorientation,  and  spatial  statistics  of  grains/phases)  on  the 
distribution  of  cyclic  and  directional  slip  and  local  stress  states  in  3D  microstructures  to  quantify  the 
probability  of  forming  and  growing  small  fatigue  cracks  [31,60-61,63-64].  The  methodology  supports  the 
notion  of  design  of  microstructures  to  enhance  fatigue  resistance  in  DOE  USCAR  [46],  DARPA  AIM 
[65-68],  and  DARPA  Prognosis  [69]  and  related  programs  [70]  that  have  utilized  this  kind  of  strategy  to 
address  sensitivity  of  fatigue  crack  formation  to  microstructures  of  engine  materials  (Ni-base  superalloys 
and  oc-|3  Ti  alloys  in  the  engine  systems  prognosis  program). 

The  primary  modeling  “workhorse”  of  these  methods  has  been  computational  crystal  plasticity,  which 
serves  to  identify  and  characterize  local  material  hot-spots  controlled  by  intrinsic  (e.g.,  local  grain 
orientations)  or  extrinsic  (e.g.,  inclusions  or  large  precipitate)  microstructure  attributes.  McDowell 
[29,60-61]  has  discussed  in-depth  the  utility  of  crystal  plasticity  in  sensitivity  studies  that  support  design 
of  microstructure  morphology  (heat  treatment,  thermo-mechanical  processing)  for  fatigue  resistance,  with 
emphasis  on  variability  in  fatigue  life  in  the  high  cycle  fatigue  regime.  The  key  elements  of  this  approach 
are  the  identification  of  key  fatigue  damage  modes  at  the  scale  of  grains  and  assignment  of  corresponding 
nonlocal  (volume  averaged)  Fatigue  Indicator  Parameters  (FIPs).  These  elements  are  employed  in  the 
estimation  of  expected  initial  crack  distributions  based  on  assumed  nucleation  and  growth  relations  at  the 
scale  of  grains,  and  limited  small  crack  growth  simulations  within  heterogeneous  microstructure. 
Different  FIPs  [1,61-62]  have  been  introduced  to  reflect  the  relative  roles  of  reversed  and  cumulative 
directional  slip  at  the  scale  of  individual  grains  or  crystalline  regions  in  polycrystals  (i.e.,  microplasticity) 
on  driving  formation  of  cracks,  either  transgranular  along  slip  bands  or  intergranular  due  to  progressive 
slip  impingement,  as  shown  in  Figure  Al.  In  high  cycle  fatigue  (HCF),  the  majority  of  component  life  is 
typically  spent  in  processes  of  forming  and  growing  cracks  in  the  microstructurally  small  crack  (MSC) 
regime,  with  crack  sizes  on  the  order  of  characteristic  microstructure  length  scales  (grain  size,  phase  or 
cluster  spacing  or  diameter).  For  example,  the  Fatemi-Socie  [71-72]  shear-based  FIP  (two  parameter 
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approach)  has  been  shown  to  correlate  multiaxial  fatigue  crack  initiation  data  very  well  in  both  LCF  and 
HCF  regimes  at  the  grain  scale  and  above  [63,73-74],  It  is  defined  by 


a  yp 

FIPfs  =  — 


1  +  K 


(1) 


V  y  j 

where  Ay^x  /  2  is  the  nonlocal  maximum  cyclic  plastic  shear  strain  averaged  over  a  finite  volume  of 


material  (on  the  order  of  several  pm3),  <7^*x  and  <jy  are  the  maximum  stress  normal  to  the  plane  of 

A/Cx  and  cyclic  yield  strength,  respectively,  and  constants  mediates  the  influence  of  normal  stress. 

The  spatial  volume  for  nonlocal  averaging  of  the  FIP  may  be  defined  according  to  the  nature  of  the 
simulation,  and  is  desirable  both  for  purposes  of  numerical  regularization  (mesh  insensitivity)  and 
targeting  length  scales  associated  with  crack  embryos  (e.g.,  slip  band  width/spacing,  inclusion  size,  etc.). 
A  similar  but  distinct  parameter  has  been  introduced  by  McDowell  [61]  to  address  correlations  of 
directional  plastic  strain  accumulation  with  crack  formation  due  to  grain  or  phase  boundary  impingement 
(Zener  pileup),  i.e., 


FIPzp=K 


C  ST*1*  ^ 

1 +  K  ^ 


v 


(2) 


y  J 


where  y?*  =ma is  the  maximum  net  plastic  shear  strain  averaged  over  a  finite  volume  of 
material  among  all  planes  with  unit  tangent  and  normal  vectors  m  and  n,  respectively,  <j^B  is  the  average 


peak  stress  normal  to  a  boundary  segment  impinged  by  this  slip,  and  K  is  a  constant  that  mediates  the 

effect  of  this  normal  stress.  This  kind  of  directional  (or  ratchet)  strain  accumulation  is  key  to 
understanding  rolling  contact  fatigue  processes  within  heterogeneous  microstructures  consistent  of  grains 
of  various  sizes,  shapes  and  orientations  and  nonmetallic  inclusions,  of  high  relevance  to  life-limiting 
situations  in  bearing  applications,  for  example. 

The  philosophy  of  this  mesoscopic  approach  to  fatigue  lifetime  prediction  and  characterization  of 
variability  of  multiple  stages  of  crack  nucleation  and  growth  has  been  elucidated  in  recent  reviews  by 
McDowell  [61-62]  and  McDowell  and  Dunne  [1].  McDowell  and  colleagues  [77]  have  developed 
schemes  by  which  microstructure-sensitive  poly  crystal  plasticity  models  for  duplex  Ti-6A1-4V  [78-81], 
y—y'  Ni-base  superalloys  [82-83],  and  martensitic  steels  [84-86]  can  be  analyzed  in  terms  of 
distributions  of  slip  and  mechanism-relevant  continuum  FIPs,  including  effects  of  pores  and  nonmetallic 
particles  on  initiation  of  fatigue  cracks.  When  employed  in  the  context  of  poly  crystal  plasticity,  these 
FIPs  enable  assessment  of  the  relative  potency  to  nucleate  and  grow  cracks  at  the  scale  of  individual 
grains  by  assigning  a  relation  between  the  number  of  cycles  necessary  to  nucleate  a  crack  at  the  scale  of 

the  nonlocal  averaging  volume  and  the  associated  FIP,  i.e.,  Nnucl  =  /  (FIP) ,  for  example  using  an 


extension  of  a  slip  band-based  nucleation  criterion  [83]. 

The  methodology  has  also  been  extended  to  a  multistage  decomposition  of  the  total  fatigue  life  given 
by  McDowell  et  al.  [46,60-62],  i.e., 


N  =  N  +  N  +  N 

iyT  ^nucl  ^  ^MSC/PSC  ^  1NLC 


(3) 


where  Nnucl  is  the  number  of  cycles  necessary  to  form  a  crack  embryo  at  scales  on  the  order  of  microns, 


N msc/psc  *s  ^e  number  of  cycles  necessary  to  grow  microstructurally/physically  small  cracks  through 

the  adjacent  microstructure,  and  NLC  is  the  number  of  cycles  of  growth  of  mechanically  long  cracks  that 

meet  similitude  requirements  of  Linear  Elastic  Fracture  Mechanics  (LEFM).  Typically,  fatigue  crack 
nucleation  occurs  crystallographically  within  a  favorably  oriented  grain/phase,  particularly  for  low  to 
moderate  stacking  fault  energy  fee  crystals  or  hep  crystals  with  a  tendency  towards  planar  slip.  A  useful 


4 


Approved  for  public  release;  distribution  unlimited 


heuristic  is  to  use  grain  level  FIPs  to  project  cycles  to  form  at  crack  at  the  scale  of  individual  grains  or 
phases,  often  tens  of  microns  in  size,  which  is  a  combination  of  nucleation  and  microstrructurally  small 
crack  (MSC)  growth  phases,  i.e.,  the  incubation  life  Ninc  =  Nnucl  +  NMSC  L  grain . 

Shenoy  et  al.  [83]  have  also  applied  a  MSC  regime  growth  law  for  Stage  I  crystallographic  growth 
that  employs  the  two  -parameter  FIP  in  Eq.  (1)  (sliding  displacement  controlled  by  cyclic  plastic  shear 
strain,  with  an  additional  effect  of  peak  normal  stress  to  this  plane),  where  the  analogy  to  mixed  mode  AJ- 
integral  [87]  and  cyclic  crack  tip  displacement  range  ACTD  is  made  (cf.  [88]),  i.e., 

=  AFS  ( Ty  F1Pfs  )  a — rijb  =  <j)  ACTD  -  r^b  (4) 

MSC 

where  AFS  is  a  constant,  Ty  is  the  critical  resolved  shear  stress  on  the  slip  plane,  and  ACTDth  is  a  threshold 

crack  tip  displacement  range  below  which  the  crack  will  not  grow.  Clearly,  relation  (4.1)  employs  a  FIP 
that  is  averaged  over  an  uncracked  domain  of  a  grain  or  phase  as  a  driving  force  parameter  when 
multiplied  by  crack  length  that  is  intended  to  parameterize  the  dependence  of  the  fine  scale  propagation 
relation  in  (4.2)  on  ACTD.  Later  discussion  will  shed  light  on  the  implied  connection  of  the  quantity 
(TyFIPFS)a  to  ACTD .  In  Eq.  (4),  an  equivalence  of  the  middle  relation  in  (4.1)  is  implied  with  regard  to 

the  relation  based  on  ACTD  via  the  right  expression  (4.2),  where  b  is  the  magnitude  of  the  Burgers  vector 
and  (|)  and  is  an  irreversibility  factor  that  controls  cyclic  crack  advance  associated  with  the  fraction  of 
process  zone  dislocations  that  do  not  return  to  the  crack  tip  during  unloading.  Parameter  §  carries  the 
“imprint”  of  atomic  scale  processes  of  new  crack  surface  creation  owing  to  slip  irreversibilities  and  can 
effects  of  interactions  with  precipitates  and  impurities.  The  threshold  r^b  characterizes  the  level  of 

irreversible  range  of  the  ACTD  below  which  the  crack  should  not  extend,  with  r\.  on  the  order  of  unity. 

This  methodology  considers  the  sequential  processes  of  crack  formation  in  hot  spot  grains  during 
cyclic  loading  and  growth  through  multiple  adjacent  grains,  with  no  attempt  to  resolve  mechanisms  in  the 
sub-micron  scale  regime  related  to  incipient  stages  of  crack  formation/nucleation.  Accordingly,  it  is 
argued  that  microstructure  at  scales  of  microns  and  above,  including  size,  shape  and  spatial  distributions 
of  grains/phases,  orientation  and  misorientation  effects,  and  extrinsic  attributes  such  as  nonmetallic 
inclusions  or  pores,  has  a  dominant  influence  on  the  variability  of  HCF  and  VHCF  fatigue  life,  which  is 
spent  dominantly  within  the  regimes  of  initial  crack  formation  and  growth  in  the  surrounding  few 
grains/phases.  For  this  reason,  the  algorithms  in  the  present  paper  do  not  consider  crack  growth  into  the 
physically  small  and  long  crack  regimes. 

Under  HCF  conditions,  the  cyclic  plastic  deformation  is  highly  heterogeneous  within  the 
microstructure;  accordingly,  this  is  the  regime  in  which  variability  and  size  effects  are  most  pronounced. 
A  strategy  for  computational  HCF  modeling  of  components  that  must  last  millions  of  cycles,  such  as 
shafts,  bearings,  and  gears,  for  example,  should  focus  on  extreme  value  statistics  of  potential  sites  for 
microplastic  strain  localization  and  fracture  that  drive  crack  formation  and  early  growth.  Moreover,  the 
issue  of  arrest  of  small  cracks  that  form  at  isolated  sites  of  cyclic  plastic  strain  intensification  with  in  the 
microstructure  is  relevant  to  fatigue  strength  thresholds.  There  are  several  types  of  thresholds  [89].  One 
is  the  absence  of  microplasticity  within  grains  sufficient  to  nucleate  cracks  or  to  drive  growth  of  micron- 
scale  embryonic  cracks  within  individual  grains/phases.  Another  somewhat  higher  threshold  is  associated 
with  crack  arrest  at  either  the  first  or  second  major  impediment  to  slip,  typically  either  grain  or  phase 
boundaries.  The  latter  is  commonly  referred  to  as  a  non-propagating  threshold  or  fatigue  limit,  and  is 
relevant  to  HCF  conditions.  The  realm  of  thresholds  in  VHCF  is  relatively  unexplored,  and  it  is 
interesting  to  consider  that  a  change  of  deformation  mechanism  or  dislocation  substructures  at  these  lower 
cyclic  plastic  strain  levels  may  give  rise  to  thresholds  for  which  we  have  little  understanding  at  present. 

Prediction  of  effects  of  variable  amplitude  load  histories  and  prediction  of  mean  fatigue  life  as  a 
function  of  microstructure  are  objectives  with  added  complexity  and  are  presently  beyond  the  purposes 
and  scope  of  the  class  of  approaches  considered  here.  In  fact,  true  prediction  of  mean  fatigue  life  requires 
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complete  understanding  of  crack  nucleation  processes,  which  involve  dislocation  substructures  and 
mechanisms  at  the  nanoscale,  including  specific  effects  of  composition,  solid  solution  strengthening,  etc. 
Fatigue  crack  nucleation  is  indeed  a  grand  challenge  problem,  since  it  considerable  idealization  may  be 
involved.  Very  recent  works  of  Sangid  et  al.  [90-91]  are  promising  in  this  regard  and  consider  not  only 
the  grain  orientation  distribution  but  also  address  the  grain  boundary  character  distribution  using  energetic 
approaches  for  nucleation  that  essentially  lend  more  physics  to  earlier  concepts  of  dislocation  dipole 
models  for  crack  nucleation  by  Tanaka  and  Mura  [92]  and  Venkarataman  et  al.  [93]  that  combine  slip 
banding  with  an  energetic  description.  They  have  also  demonstrated  the  ability  to  project  variability  of 
HCF  life  for  a  polycrystalline  ensemble.  We  note,  however,  that  there  are  yet  unresolved  questions 
regarding  the  role  of  vacancies  in  the  emergent  instability  of  a  crack  forming  in  a  lattice.  Work  in  the  role 
of  generation  of  point  defects  (mainly  vacancies)  due  to  production  and  annihilation  of  dislocations  within 
regions  of  localized  cyclic  plastic  deformation  goes  back  several  decades  [94-96].  Using  positron 
annihilation,  Egger  et  al.  [97-98]  have  shown  that  vacancy  concentration  evolves  in  fatigue  (but  not 
monotonic  fracture)  within  regions  of  localized  slip  to  the  point  at  which  cracks  nucleate.  The 
significance  of  these  studies  is  that  the  driving  force  for  crack  nucleation  is  affected  by  the  production  and 
migration  of  point  defects  and  not  just  based  on  dislocation  pileups,  as  proposed  by  Tanaka  and  Mura 
[92]. 

For  the  time  being,  we  proceed  in  microstructure-sensitive  design  of  advanced  alloys  without  explicit 
predictive  models  for  nucleation,  owing  to  the  role  of  microstructure  in  affecting  small  fatigue  crack 
growth  and  inducing  fatigue  limits  through  blockage  of  small  cracks.  This  can  be  augmented  in  time 
gaining  a  more  predictive  footing  with  regard  to  crack  nucleation  processes.  Of  course,  an  interesting 
question  is  what  minimal  set  of  microstructure  information  is  necessary  to  predict  variability  of  fatigue 
life  for  a  given  microstructure  and  loading  condition.  It  is  not  fully  resolved,  but  suffice  it  to  say  that 
attention  is  drawn  to  the  life  limiting  features  in  HCF  such  as  largest  favorably  oriented  grains,  adjacent 
grains  with  favorable  orientation,  interactions  among  grains  that  enhance  slip,  nonmetallic  inclusions  at 
higher  scales,  etc.  Mesoscale  simulations  of  the  type  considered  here  can  provide  much  useful 
information  even  in  the  absence  of  predictive  nucleation  relations  that  reflect  composition  and  grain 
boundary  character. 

Quantifying  the  sensitivity  of  fatigue  life  to  microstructure  opens  several  avenues: 

•  The  ability  to  tailor  microstructure  to  improve  fatigue  resistance 

•  Comparison  of  candidate  microstructures  for  specific  applications,  including  low  probability 
of  failure  regimes  of  interest  to  design  applications 

•  The  ability  to  augment  costly  experiments  in  the  HCF  and  VHCF  regimes  to  quantify 
variability  and  scatter  in  fatigue 

This  paper  will  consider  mainly  the  latter  two  items  in  the  list  above.  However,  they  enable  the  first 
one,  namely  providing  decision-support  for  tailoring  microstructure  to  achieve  enhanced  fatigue 
resistance.  Figure  2  shows  a  loop  in  microstructure  design  that  employs  the  microstructure-sensitive 
simulation  methods  outlined  here.  By  subjecting  digital  statistical  volume  elements  (SVEs)  constructed 
in  adherence  with  spatial  statistics  of  actual  microstructures  to  computational  simulation  for 
representative  loading  histories,  a  distribution  of  microstructure  (e.g.,  grain)  scale  FIPs  are  computed  for 
many  SVEs  and  the  histogram  used  to  assess  the  range  of  expected  minimum  fatigue  lives  for  a  given 
nominal  microstructure  and  loading  condition.  This  information  is  then  fed  back  into  the  instantiation  of 
digital  SVEs  with  modified  microstructure  aimed  at  increasing  minimum  life  or  reducing  variability, 
according  to  design  targets.  There  are  several  potential  design  applications  of  this  kind  of  approach.  For 
example,  (i)  designing  application-specific  composition  and  process  route,  (i)  grading  microstructures  in 
components  (location  specific  properties),  and  (iii)  accelerating  certification  of  design  for  minimum  life. 

Within  this  paradigm,  there  are  several  important  points.  First,  it  is  assumed  that  variability  emerges 
mainly  from  the  state  of  cyclic  inelastic  deformation  within  individual  grains/phasees,  which  is  affected 
primarily  by  grain  size,  shape,  orientation  and  relation  to/coupling  with  neighboring  grains.  This  is 
particularly  important  in  HCF  since  only  favorable  grains  are  activated  within  a  highly  heterogeneous 
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field.  Second,  the  constitutive  models  for  deformation  and  for  relating  FIPs  to  lifetime  estimates  should 
be  calibrated  to  mean  fatigue  response  of  specific  microstructures  (validation),  for  several  reasons: 

•  uncertainty  in  idealization  and  approximation  of  models  (e.g.,  slip,  nucleation,  crack  growth) 

•  potential  omission  of  deformation  or  fatigue  damage  mechanisms 

It  may  seem  counter-intuitive  but  absolute  prediction  of  fatigue  life  is  well  beyond  current  capabilities 
of  this  kind  of  modeling  framework  until  relations  for  fatigue  crack  nucleation  and  microstructurally 
small  crack  growth  are  better  understood  and  quantified.  When  mean  fatigue  response  is  used  to  calibrate 
models  for  a  given  fatigue  loading  condition  and  set  of  microstructures,  so  long  as  the  dominant 
mechanisms  don’t  change,  it  is  assumed  that  mean  fatigue  responses  can  be  projected  for  the  same 
loading  conditions  and  other  microstructures  within  the  range  calibrated  or  modeled. 

3.  Multiaxial  Fatigue  and  FIPs 

One  of  the  primary  motivating  features  for  using  FIPs  is  their  simplicity  as  computable  mesoscopic 
parameters,  while  demonstrating  desirable  predictive  character  for  multiaxial  stress/strain  states  of  the 
kind  that  grains/phases  actually  experience.  A  wide-range  of  FIPs  have  been  proposed  in  the  last  30  years 
to  characterise  the  formation  and  early  growth  of  fatigue  cracks.  Socie  and  Marquis  [100]  summarized 
many  of  the  FIPs  devised  on  stress,  strain  or  energy  concepts,  which  supports  that  a  single  continuum 
parameter  cannot  correlate  equally  well  with  the  multiaxial  fatigue  crack  initiation  behavior  of  all  metals. 
For  example,  the  FIP  proposed  by  Fatemi  and  Socie  [71]  in  Eq.  (1)  based  on  the  critical  plane  approach 
correlates  well  with  the  small  fatigue  crack  formation  and  growth  in  ductile  metals  that  sustain  shear- 
based  fatigue  damage. 

3.1  Coarse  Scale  (Uncracked)  FIP  Correlations  for  Multiaxial  Fatigue 

As  suggested  by  Reddy  et  al.  [101]  and  McDowell  and  Berard  [88],  the  Fatemi-Socie  FIP  can  be 
understood  as  a  local  fatigue  driving  force  similar  to  the  AK  or  the  AJ  customized  to  describe  the  early 
crystallographic  crack  growth  phenomena.  Indeed,  McDowell  and  Berard  [88]  and  McDowell  [63] 
showed  that  the  Fatemi-Socie  parameter  in  Eq.  (1)  resulted  in  fatigue  life  correlations  over  a  broad  range 
of  multiaxial  loading  conditions  in  a  manner  similar  to  small  fatigue  crack  growth  relations  motivated  by 
the  AJ  -integral  of  Elastic-Plastic  Fracture  Mechanics  (EPFM).  Hoshide  and  Socie  [87]  appealed  to  the 
elastic  and  plastic  forms  of  the  standard  J- integral  expression  of  EPFM  to  model  mixed  mode  crack 
growth  in  multiaxial  fatigue,  with  J  given  by 

7r(<J2+T2) 

J  —  —  aejj  +  /  ( n,  )  G£  a  —  J e  +  J p  (5) 

Here,  a  and  aeff  are  actual  and  effective  crack  lengths,  and  the  stress  biaxiality  ratios  are  given  by 

X0—T I  <Jyy  and  =  <Jxx  /  Gyy ,  where  T  and  Gyy  are  the  far  field  shear  and  crack  plane  normal  stresses, 

respectively,  and  (Jxx  and  is  the  direct  stress  parallel  to  the  crack.  Considering  AJ  for  cyclic  loading 

based  on  the  loading  part  of  the  cycle,  they  correlated  crack  growth  data  for  IN  718  for  cracks  less  than  1 
mm  in  length  using  the  power  law  relation 


da  /  A 

- =  Cj(AJ)  (6) 

dN  JK  ’ 

where  C3  and  M  3  depend  on  the  biaxiality  ratios,  and M  3  ranged  from  1.31  to  1.45.  It  turns  out  that 
small  crack  growth  data  (crack  lengths  on  order  of  hundreds  of  microns,  beyond  the  follow  reasonably 
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well  this  kind  of  correlation  based  on  the  homogeneous,  macroscopic  value  of  AJ  (cf.  McDowell, 
1996HCF).  This  is  perhaps  not  surprising  owning  to  the  linear  relation  between  J  and  CTOD  even 
under  EPFM  conditions  [102]. 

By  effectively  asserting  a  coarse  scale  (based  on  uncracked  cyclic  plastic  strain  field)  FIPfs  -based 

growth  law  of  the  type  given  in  Eq.  (4)  for  physically  small  cracks  under  plasticity-dominated  LCF 
conditions,  McDowell  [63]  showed  that  the  characteristic  contours  for  multiaxial  fatigue  in  the  so-called 

rp  -plane  were  in  qualitative  agreement  with  trends  from  experimental  data,  as  shown  in  Fig.  3.  TheT  - 
plane  was  introduced  by  Brown  and  Miller  [103]  and  presents  trajectories  of  constant  fatigue  life  data 
corresponding  to  a  crack  of  a  given  length,  say  500  pm,  in  terms  of  maximum  range  of  shear  strain  versus 
peak  normal  strain  to  the  same  plane.  McDowell  and  Berard  [88]  introduced  a  more  complex  variant  of 
this  growth  relation  with  more  degrees  of  freedom  to  match  experimental  data  for  each  of  the  two 
branches  evident  in  Fig.  3. 

Alternatively,  when  the  nonlocal  FIP  is  computed  in  the  vicinity  of  the  tip  of  a  small  crack,  it  includes 
the  plastic  strain  intensification  due  to  crack  and  the  relation  in  Eq.  (4)  may  written  as 

=  AFSFIPFS,crackd-ACTDth  (7) 

MSC 

where  d  is  a  characteristic  length  scale  employed  for  averaging  the  fine  scale  FIP.  In  this  case,  the 
product  FIPfs  crack  d  is  intended  to  parameterize  theACTD;  for  this  to  hold,  it  must  turn  out  that 

FIPfs  crack  is  linearly  related  to  ACTD .  We  will  explore  the  efficacy  of  FIPfs  crack  in  Eq.  (7)  in  the  next 
section,  with  reference  to  driving  force  for  crystallographic  cracks  in  single  crystals. 

3.2  Microstructurally  Small  Crack  Growth  -  Direct  Relation  of  FIPfs,  crack  to  ACTD 

In  this  section  we  compare  the  Fatemi-Socie  FIP  with  the  damage-tolerant  approaches  (e.g.,  LEFM, 
EPFM)  have  traditionally  estimated  the  fatigue  crack  growth  driving  force  by  means  of  the  range  of  stress 
intensity  factor  (AK),  the  range  of  J-integral  (AJ)  or  the  range  of  crack  tip  displacement  (ACTD).  The 
experimental  validation  of  these  magnitudes  is  usually  based  on  far-field  measurements  (an  exception  is 
Schwalbe’s  b5  [104])  that  rely  on  the  similitude  assumption  of  fracture  mechanics.  However,  as  discussed 
by  McDowell  [29],  the  similitude  assumption  is  often  not  valid  when  the  crack  length  or  cyclic  plastic 
zone  size  is  on  the  order  of  the  characteristic  dimensions  of  the  microstructure.  In  this  case, 
microstructure-sensitive  characterizations  need  to  consider  local  fields  to  estimate  fatigue  driving  forces. 
Since  the  concept  of  a  ACTD  does  not  rely  on  the  similitude  hypothesis,  it  is  a  potential  measure  of  the 
local  crack  driving  force,  but  requires  detailed  FEM  simulations  to  capture  the  influence  of  load  sequence 
and  microstucture  in  the  vicinity  of  the  crack  tip.  Recently,  Castelluccio  and  McDowell  [105-106] 
compared  the  fine  scale  Fatemi-Socie  FIPfs  crack  (i.e.,  computed  within  a  cracked  single  grain  or  crystal) 

to  the  ACTD  using  finite  element  simulations  of  single  crystals  with  cracks  on  {111}  slip  planes  under 
shear  and  mixed  mode  loading.  For  simplicity,  we  omit  the  subscript  for  the  remainder  of  this  section  and 
refer  to  FIPfs  crack  simply  as  “FIP“.  They  employed  3D  crystal  plasticity  constitutive  models  within  user 

material  subroutines  for  the  FE  code  ABAQUS  [107]  for  both  OFHC  Cu  and  RR1000  Ni-base  superalloy 
single  crystals  oriented  for  single  slip  in  remote  shear  loading,  with  a  sense  parallel  to  the  crack.  A  total  of 
five  models  with  stationary  crack  lengths  of  2  pm,  5  pm,  7.5  pm,  10  pm  and  15  pm  were  analyzed;  Fig.  4 
shows  one  example.  The  upper  and  lower  faces  of  the  crystal  in  Fig.  4  were  cyclically  displaced  in  shear 
or  mixed  mode  loading,  as  shown  in  Fig.  5;  three  loading  cycles  were  applied  with  the  magnitude  of  the 
displacement  vector  D  ranging  from  0.005  pm  to  0.1  pm,  equivalent  for  this  geometry  to  nominal  peak 
strains  from  0.05%  to  1%.  The  strain  ratio  R£  =  £min/£max  —  0  was  employed,  where  £max  and  £min  are 
the  maximum  and  minimum  applied  strain  over  a  loading  cycle.  For  homogeneous  single  crystals,  these 
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values  resulted  in  plastic  strain  amplitudes  that  spanned  the  HCF  to  LCF  regimes.  The  FIP  was  averaged 
in  this  case  within  a  band  of  uniform  2  pm  thickness  parallel  to  the  crack,  as  shown  to  the  right  in  Fig.  4. 

The  ACTD  was  evaluated  by  measuring  the  displacement  between  red  nodes  in  Fig.  4  at  the  mid¬ 
thickness  over  the  loading  portion  of  the  third  loading  cycle.  The  total  crack  tip  displacement  range  was 
calculated  as 


ACTD= V  ACTSD2  +ACTOD2  (8) 

based  on  the  work  of  Ma  et  al.  [108]  and  employed  by  Bennett  and  McDowell  [39-41].  The  cyclic  crack 
tip  opening  (ACTOD)  and  sliding  (ACTSD)  displacements  were  computed  based  on  the  change  in 
distance  between  red  nodes  shown  in  Fig.  4. 

Because  the  ACTD  includes  elastic  and  plastic  contributions,  the  FIP  was  calculated  for  each  element 
based  on  the  total  strain  tensor.  Figure  6  compares  the  ACTD  with  the  fine  scale  FIP  for  shear  (left)  and 
mixed  mode  (right)  loading  and  OFHC  Cu  and  RR1000  Ni-  base  superalloy.  Even  though  the  averaging 
volume  was  defined  in  terms  of  a  band  as  shown  in  Fig.  4,  Cu  and  RR1000  exhibit  similar  behavior  with 
regard  to  the  correlation  of  the  fine  scale  FIP  with  ACTD  under  single  slip  shear  loading.  The  correlation 
for  Cu  under  mixed  mode  loading  is  nearly  identical,  whereas  for  the  Ni-base  superalloy  there  is 
increased  variability  for  mixed  mode  loading.  Results  in  Fig.  6  essentially  demonstrate  linearity  of  the 
ACTD  with  the  FIP  for  both  loading  modes  for  ACTD  values  above  1  nm,  which  is  on  the  order  of  the 
magnitude  of  the  Burgers  vector. 

Hence,  it  is  interesting  to  note  that  whether  a  coarse  scale  (uncracked)  or  fine  scale  (cracked)  FIPfs  is 
employed,  there  is  close  respective  connection  to  both  multiaxial  fatigue  behavior  of  shear-dominated 
crack  initiation  (formation  and  early  growth)  processes  at  the  mesoscale  [63,88]  and  the  ACTD  driving 
force  for  crystallographic  fatigue  crack  growth  in  single  crystals.  This  suggests  utility  of  the  FIP  concept 
from  single  grains  to  mesoscales  involving  many  grains,  regardless  of  whether  similitude  requirements  of 
LEFM  are  met. 

4.  Microstructure  Attributes  and  Extreme  Value  Correlations  in  HCF 

4.1.  Grains  (size,  shape,  orientation,  and  disorientation  distributions) 

Grain  size  can  have  a  significant  effect  on  the  fatigue  failure  mechanism  in  Ti  alloys  and  Ni-base 
superalloys  with  a  given  microstructure.  The  grain  boundary  network  is  also  relevant.  Directionally 
solidified  and  single-crystal  forms  of  Ni-base  superalloys  suppress  or  eliminate  grain  boundary 
diffusional  creep  and  intergranular  oxidation  mechanism  at  elevated  temperatures.  The  AIM  program 
[65-67]  hinted  at  the  possibility  of  location  specific  design  of  components  (e.g.,  dual  or  hybrid 
microstructures)  with  optimized  microstructure  to  meet  operational  requirements  (stress  state, 
temperature,  environment,  etc.)  in  different  regions  of  components. 

The  mode  of  fatigue  crack  formation  and  growth  in  Ni-base  superalloys  at  elevated  temperatures 
strongly  depends  on  the  grain  size  of  the  microstructure.  For  subsolvus,  fine  grain  (average  grain  sizes  - 
3-15 fim)  Ni-base  superalloys,  cracks  tend  to  form  via  a  Stage  II,  non-crystallographic,  transgranular 
cleavage  mechanism  near  non-metallic  inclusions  (NMIs)  or  pores  [109-111];  this  fatigue  crack  growth 
mechanism  transforms  to  intergranular  growth  for  surface  cracks  when  the  crack  reaches  a  characteristic 
length  (-50-100 jum)  [109-111]  or  when  dwell  periods  occur  [110];  for  sub-surface  cracks,  intergranular 
cracking  requires  crack  lengths  greater  than  the  characteristic  length  and  exposure  of  the  crack  to  the 
surface  [110].  Alternatively,  due  to  the  planar  slip  mechanism  in  Ni-base  superalloys  [111-112], 
supersolvus,  coarse  grain  (average  grain  sizes  >  15 /urn )  Ni-base  superalloys  favor  Stage  I  crystallographic 
crack  formation  from  inclusions,  pores  or  large  grains.  This  fatigue  mechanism  has  been  observed  in 
Rene’  88DT  [113-115],  IN100  [116],  Waspoloy  [109],  and  HIP  Astroloy  [109].  The  prediction  of  the 
interplay  between  (1)  the  non-crystallographic,  transgranular  cleavage  failure  mechanism  near  NMIs  or 
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pores  for  fine  grain  microstructures  and  (2)  the  crystallographic  crack  formation  failure  mechanism  in 
coarse  grain  microstructures  becomes  very  important  in  the  design  of  dual-microstructure  Ni-base 
superalloy  disks  (e.g.,  [117]). 

Another  mesoscopic  parameter  that  affects  fatigue  significantly  is  the  distribution  of  grain  orientation, 
i.e.,  the  texture  of  the  material.  Texture  development  becomes  especially  important  for  highly  anisotropic 
materials  such  as  the  HCP  Ti  alloys  modeled  in  this  study.  For  example,  in  duplex  a  I  /?Ti-6Al-4V  alloys 
slip  is  preferred  in  the  primary  a  phase  via  basal  and  prismatic  slip  planes  since  these  slip  planes  have  a 
lower  critical  resolved  shear  stress  as  compared  to  the  primary  and  secondary  pyramidal  slip  planes  [78- 
81,118].  The  fact  that  slip  is  easier  on  basal  and  prismatic  planes  can  have  major  implications  for  rolled 
Titanium  microstructures  or  clusters  of  grains  that  are  preferentially  oriented  for  basal  or  prismatic  slip. 

The  minimum  misorientation  (i.e.,  disorientation)  between  neighboring  grains  plays  an  important  role 
in  fatigue  crack  formation  and  growth.  Misorientation  is  the  relative  change  in  orientation  between 
neighboring  grains  and  can  be  denoted  by  a  rotation  angle/axis  pair  [119-120].  Grain  boundaries  with  a 
high  disorientation  angle  are  often  more  favorable  locations  for  fatigue  crack  formation  [121].  This  is  due 
to  the  progressive  pile-up  of  dislocations  in  slip  bands  (Zener  mechanism)  during  cyclic  loading  that 
impinge  on  grain  or  phase  boundaries,  or  upon  oxidized  inclusion  interfaces,  which  can  lead  to  formation 
and  propagation  of  small  cracks  in  the  microstructure.  Alternatively,  high  angle  grain  boundaries  tend  to 
hinder  or  arrest  crack  propagation.  Specifically,  the  tilt  and  twist  angles  of  a  grain  boundary  directly 
affect  the  ability  to  transmit  plasticity/slip  across  the  grain  boundary.  Moreover,  twin  boundaries  play  a 
significant  role  in  fee  systems  [122]. 

Experimental  techniques  combining  various  technologies  such  as  serial  sectioning,  X-ray  diffraction 
contrast  tomography,  EBSD/OIM,  and  digital  image  correlation  allow  for  very  accurate  characterization 
of  microstructure  statistical  information  including  grain  shape,  grain  size,  grain  orientation  and  grain 
misorientation  [123-125].  Although  collection  of  this  high  resolution  stochastic  information  is  crucial  for 
more  accurate  reconstruction  and  modeling  of  realistic  3D  microstructures,  these  techniques  can  be 
extremely  time-consuming.  This  makes  it  experimentally  infeasible  to  obtain  this  level  of  detail  about 
every  simlation  that  we  run.  Thus,  we  use  the  statistical  information  obtained  from  limited  number  of 
experimental  characterizations  as  input  parameters  for  the  computational  reconstruction  and  simulation  of 
statistically  representative  3D  microstructures. 

We  employ  two  different  means  of  polycrystalline  reconstruction,  a  traditional  Voronoi  tessellation 
technique  and  a  random  sequential  adsorption  (ellipsoidal  packing)  algorithm.  The  Voronoi  tessellation 
method  is  used  to  create  microstructure  instantiations  with  convex  grain  hulls.  The  fitting  of  these 
microstructures  to  a  targeted  grain  size  distribution  often  requires  a  secondary  Monte  Carlo  or  simulated 
annealing  algorithm  to  optimize  the  grain  size  distribution.  The  use  of  Voronoi  tessellation  method  has 
some  limitations.  First  of  all,  there  is  a  limited  ability  to  be  able  to  reproduce  bimodal  or  multimodal 
grain  size  distributions  126].  In  addition,  the  average  number  of  neighboring  cells  in  a  random  Voronoi 
tessellation  can  differ  from  that  seen  in  experiments.  For  example,  Groeber  et  al.  [124-125]  showed  that 
the  average  number  of  neighboring  cells  in  a  fine-grained  IN100  Ni-base  superalloy  was  12.9,  which 
differs  from  that  predicted  (15.5)  for  a  random  Voronoi  tessellation  [126].  One  method  to  overcome  these 
drawbacks  of  the  Voronoi  tessellation  method  is  to  use  a  random  sequential  adsorption  method  using 
ellipsoids  to  represent  grains.  The  approximation  of  grains  as  ellipsoids  has  been  undertaken  previously 
by  many  researchers  [123-124,127-129].  In  Refs.  [127-128],  EBSD  analysis  of  orthogonal  serial  sections 
was  used  to  create  ellipsoidal  grains  and  reconstruct  the  polycrystalline  microstructure  for  a  hot  rolled 
pure  aluminum.  Similarly,  Groeber  et  al.  [123-124]  used  serial  sectioning  to  characterize  and  reconstruct 
statistically  representative  Ni-base  superalloy  IN  100  microstructures  that  incorporate  distributions  of 
grain  size,  grain  orientation  and  grain  misorientation.  Przybyla  and  McDowell  [129]  used  grain 
equivalent  ellipsoids  to  more  accurately  model  the  bimodal  grain  distribution  of  a  duplex  Ti-6A1-4V 
microstructure. 
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4.2.  Nonmetallic  Inclusions  (NMIs)  and  Pores 

Inclusions  and  pores  are  introduced  into  powder  metallurgy  (PM)  Ni-base  superalloys  during 
materials  processing.  As  PM  processing  techniques  improve,  cleaner  Ni-base  superalloys  are  being 
developed  with  reduced  number  densities  of  inclusions.  However,  proper  probabilistic  modeling  of  PM 
Ni-base  superalloy  components  should  take  into  account  the  probability  of  occurrence  of  a  life-limiting 
NMI  and  the  effect  of  the  inclusion  on  the  overall  probability  of  failure.  This  approach  should  also 
consider  the  size  effects  associated  with  testing  specimen-sized  versus  component-sized  volume  domains, 
since  the  statistically  weakest vv defect"  size  within  a  volume  can  change  with  the  size  of  the  volume  [ISO- 
131].  In  the  absence  of  a  NMI,  cracks  can  be  initiated  at  pores  located  near  the  surface;  Porter  III  et  al. 

[132]  performed  load  control  fatigue  testing  on  two  PM  supersolvus  Ni-base  superalloys,  Rene'  88DT  and 
IN100,  with  average  grain  sizes  of  30  and  25  pm,  respectively.  For  Rene'  88DT,  cracks  initiated 
primarily  at  non-metallic  inclusions  (NMI)  near  the  specimen  surface.  For  IN100,  life-limiting  cracks 
initiated  from  pores  70%  of  the  time  and  NMI  the  other  30%.  Also,  the  cracks  initiated  most  frequently 
at  the  surface  of  the  specimen  (60%). 

The  size  of  the  inclusion  relative  to  the  statistically  largest  grain  size  requires  further  explanation. 
Gabb  et  al.  [117]  investigated  the  fatigue  mechanisms  of  a  graded,  dual-microstructure  LSHR  Ni-base 
superalloy  used  for  a  dual  microstructure  disk.  They  tested  specimens  from  the  inner  ring  (fine-grained), 
outer  ring  (coarse-grained),  and  the  transition  zone  (contained  both  fine  and  coarse  grains).  They  found 
that  fine  grain  inner  ring  specimens  predominantly  failed  due  to  internal  NMIs  that  were  larger  than  the 
average  grain  size,  whereas  the  coarse  grain  outer  ring  and  the  transition  zone  failed  from  coarse  grain 
facets.  Moreover,  if  the  inclusion  is  smaller  than  the  grain  size,  it  seems  to  have  little  effect  on  the 
amplification  of  driving  force  for  fatigue  crack  initiation  and  early  growth  as  compared  to  the  statistically 
largest  grain  in  the  microstructure.  On  the  other  hand,  inclusions  larger  than  the  average  grain  size  will 
most  likely  serve  as  the  fatigue-critical  hotspots  in  the  microstructure.  Thus,  the  size  of  the  largest 
inclusion  as  compared  to  the  largest  grain  size,  often  referred  to  as  the  ALA  (as  large  as)  grain  size,  is  a 
very  important  attribute  to  monitor  and  control  during  processing.  This  is  the  reason  why  control  of 
inclusion  (and  pore)  size  is  crucial  in  fine  grained  Ni-base  superalloy  microstructures  which  have  average 
grain  sizes  ~  3-5  pm;  cracks  tend  to  form  in  these  fine  grained  microstuctures  via  a  Stage  II,  non- 
crystallographic,  transgranular  cleavage-like  mechanism  near  non-metallic  inclusions  (NMIs)  or  pores 
[109-111]. 

Another  important  aspect  is  the  inclusion/matrix  interface.  We  consider  three  different  conditions: 
the  inclusion  is  completely  intact  with  the  matrix,  the  inclusion  is  cracked,  or  the  inclusion  is  partially 
debonded  with  the  matrix.  Previous  work  on  martensitic  gear  steels  showed  that  intact  inclusions  show 
negligible  driving  force  for  crack  nucleation  in  the  HCF  regime  [86].  Prasannavenkatesan  et  al.  [86] 
showed  that  partially  debonded  inclusions  provide  the  worst  case  condition  for  fatigue  crack  initiation  and 
propagation,  which  is  consistent  with  other  investigations  in  a  directionally  solidified  Ni-base  superalloy 

[133]  and  a  cast  A356-T6  aluminum  alloy  [46,134].  Therefore  we  consider  that  modeling  inclusions  as 
initially  partially  debonded  is  suitable  for  minimum  life  designs.  However,  this  interface  assumption  may 
be  overly  conservative  depending  on  the  loading  condition,  location  of  the  inclusion  (within  the 
subsurface),  and  surface  treatment. 

4.3  Extreme  Value  Marked  Correlation  Functions 

Like  most  damage  processes,  fatigue  is  often  highly  dependent  on  the  extreme  value  distributions  of 
the  microstructure  attributes  most  relevant  to  the  operative  damage  mechanisms.  For  example,  it  has  been 
shown  that  in  ductile  polycrystalline  metals  the  processes  of  crack  formation  and  propagation  are 
dominated  by  the  largest  non-metallic  inclusions  [135-137]  and  largest  grains  [115,  138-139].  Moreover, 
multivariate  extreme  value  distributions  are  likely  important  when  multiple  microstructure  attributes 
interact  to  influence  fatigue  damage  formation.  For  example,  multiple  interacting  grains  or  grain  clusters 
[140-143]  and  large  textured  macrozones  (i.e.,  microtextures)  [75,144-145]  have  been  observed  to 
influence  the  fatigue  life  in  Ti  alloys.  The  extreme  value  marked  correlation  function  (EVMCF)  [99,146] 
offers  a  quantitative  measure  to  link  the  influence  of  the  distributions  of  multiple  interacting 
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microstructure  attributes  to  the  extreme  value  response  based  on  a  given  applied  loading  condition  (e.g., 
strain  state  and  strain  amplitude).  In  practice,  the  EVMCF  is  constructed  by  quantifying  the  correlated 
microstructure  attributes  that  are  relevant  to  the  operant  mechanism  of  damage  formation  and  that  are 
coincident  with  (i.e.,  marked  by)  the  regions  associated  with  the  pertinent  extreme  value  response 
parameters,  for  example  the  FIP  in  Eq.  (1).  These  extreme  value  marked  correlations  can  be  compared  to 
the  unmarked  correlation  functions  of  the  same  type  that  are  sampled  over  the  entire  microstructure 
ensemble  to  determine  the  particular  correlated  microstructure  attributes  that  are  most  significant  relative 
to  the  extreme  value  response.  In  general,  this  methodology  can  be  applied  to  quantify  the  correlated 
microstructure  attributes  that  are  most  likely  to  be  associated  with  the  extreme  value  distributions  of 
fatigue  response,  providing  feedback  to  materials  development  as  shown  in  Fig.  7  and  described  in  the 
following. 

To  define  the  EVMCF  it  is  first  helpful  to  define  the  probability  distribution  of  an  arbitrary  extreme 
valued  response  parameter  a ,  which  can  be  defined  as  /max(a)(al^)  f°r  a  given  window  fl.  In  words, 
/max(a)  describes  the  probability  that  the  response  parameter  a  is  maximum  for  a  given  window 

(e.g.,  area,  volume)  fl.  The  EVMCF  then  can  be  defined  generally  for  n  arbitrary  correlated 
microstructure  attributes  (/J1,/?2,  ...,/?n)  as  ,  fi2 ,  ...,/?n|r1,r2,  ...rn,  H),  which  describes  the 

probability  that  n  correlated  attributed  (/J1,/?2,  ...,/?n)  separated  by  the  vectors  (r1  ,r2 ,  ...rn)  can  be 
found  in  the  same  neighborhood  as  the  extreme  value  response  a  for  a  given  window  fl.  In  practice,  two 
correlated  microstructure  attributes  (/71,/?2)  can  be  described  by  several  different  standard  correlation 
functions  such  as  the  2-point  correlation  function,  the  radial  (i.e.,  pair)  correlation  function,  the  lineal  path 
correlation  function,  nearest  neighbor  correlation  function,  etc.  Higher  correlation  functions  can  also  be 
envisioned  (e.g.,  3-point)  to  explore  correlations  between  more  than  two  microstructure  attributes.  The 
type  of  correlation  function  that  should  be  used  depends  on  the  types  of  microstructure  relationships  that 
are  most  pertinent  to  the  operant  damage  mechanisms  of  interest.  For  example,  the  radial  correlation 
function  would  not  describe  the  anisotropy  that  could  be  described  using  the  2-point  correlation  function. 

Two  targeted  studies  were  performed  recently  to  quantify  the  influence  of  poly  crystalline 
microstructure  on  fatigue  of  Ni-base  superalloy  IN100  [141]  and  Ti-6A1-4V  [142,147]  using  the  EVMCF. 
The  microstructure  attribute  considered  was  the  apparent  Schmid  factor  (calculated  based  on  the 
orientation  of  the  slip  plane  relative  to  the  loading  direction).  The  response  was  defined  in  terms  of  the 
driving  forces  for  fatigue  crack  formation  as  predicted  by  the  Fatemi-Socie  FIP  in  Eq.  (1).  Multiple  SVEs 
were  instantiated  based  on  the  microstructure  statistics  of  experimentally  characterized  microstructure 
volumes  and  simulated  using  physics  based  crystal  plasticity  models  calibrated  to  experiments.  The  SVEs 
were  meshed  and  simulated  using  the  FEM  with  periodic  boundary  conditions  to  simulate  subsurface 
fatigue  crack  formation  and  in  some  cases  traction  free  boundary  conditions  on  some  surfaces  to  simulate 
surface  fatigue  crack  formation.  Correlations  between  the  apparent  Schmid  factors  were  evaluated  in 
terms  of  the  radial  correlation  function  such  that  Rmax^(jn9  ,m9'\r,  £l)dr  is  the  probability  that  the 
apparent  Schmid  factor  m9  for  grain  orientation  g  is  located  coincident  with  the  extreme  value  FIP  P  in 
H  and  within  a  distance  of  r  to  r  +  dr  of  the  apparent  Schmid  factor  m9'  associated  with  grain 
orientation  g'  in  any  specified  direction.  Calculating  the  extreme  value  distribution  of  the  FIPs  in 
conjunction  with  the  extreme  value  marked  correlation  function  provides  a  description  of  the  extreme 
value  response  of  the  microstructure  as  represented  by  the  FIP  P  and  considers  the  probabilities  of 
correlated  Schmid  factors  m9  and  m9'  (associated  with  the  grain  orientations  g  and  g')  existing  in  the 
proximity  of  extreme  values  of  P  in  a  microstructure  window  H.  By  comparing  Rm2LX(p)(jn9 ,m9'\r,  £l)dr 
with  the  correlation  function  for  the  same  Schmid  factors  (i.e.,  grain  orientation)  sampled  from  the 
complete,  bulk  material  ensemble,  given  by  R(m9 ,m9'\r,  H)dr,  the  Schmid  factors  that  are  most 
probable  to  exist  at  locations  of  extreme  value  response  can  be  identified. 

Figure  8  shows  the  extreme  value  distributions  of  the  grain  scale-averaged  Fatemi-Socie  FIP  (Eq.  (1)) 
for  three  different  magnitudes  of  the  applied  maximum  cyclic  plastic  strain  after  three  cycles  for 
simulations  of  the  PM  Ni-base  superalloy  IN100  [141].  In  these  simulations,  periodic  boundary 
conditions  were  applied  in  all  directions  to  simulate  a  bulk  subsurface  material  condition  under  loading  in 
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strain  control  with  completely  reversed  strain-controlled  loading  (i.e.,  R£  =  — 1).  As  can  be  seen  in  Fig.  8, 
in  all  cases  the  extreme  value  distributions  of  the  Fatemi-Socie  FIPs  were  well  described  by  the  Gumbel 
extreme  value  distribution  (i.e.,  R2  >  0.97).  It  is  noted  here  that  these  results  do  not  consider  the  presence 
of  voids  or  non-metallic  inclusions,  even  though  it  is  recognized  that  non-metallic  inclusions  play  a  very 
important  role  in  the  formation  of  life-limiting  fatigue  cracks  in  this  particular  material  system.  However, 
the  understanding  of  extreme  value  slip  behavior  is  very  important  in  contributing  to  the  understanding  of 
fatigue  crack  formation  in  the  vicinity  of  inclusions/voids. 

Figure  9  shows  the  extreme  value  marked  radial  distributions  for  the  apparent  Schmid  factor 
associated  with  SVEs  subjected  to  an  applied  strain  amplitude  of  0.005.  Both  octahedral  (110){111}  and 
cube  slip  (110){100)  were  modeled  explicitly  in  the  crystal  plasticity  formation  and  high  Schmid  factors 
for  each  were  compared.  Here  we  observed  that  grains  oriented  unfavorably  for  octahedral  slip  or 
favorably  oriented  for  cube  slip  are  predicted  to  exist  with  high  probability  at  the  locations  of  extreme 
value  fatigue  response  in  IN  100.  Moreover,  clusters  of  grains  oriented  for  cube  slip  or  clusters  of  grains 
oriented  for  cube  slip  surrounded  by  grains  oriented  favorably  for  octahedral  slip  are  predicted  to  exist 
with  high  probability  at  the  location  of  the  extreme  value  fatigue  response.  These  simulations  support  the 
observations  made  in  IN100  by  Li  et  al.  [116]  who  observed  fatigue  damage  formation  along  {100} 
planes  in  grains  oriented  unfavorable  for  octahedral  slip.  This  also  supports  the  hypothesis  that  the  driving 
force  for  fatigue  crack  formation  is  highest  in  clusters  of  grains  oriented  for  cube  slip  near  grains  oriented 
for  octahedral  slip.  In  the  presence  of  non-metallic  inclusions,  it  would  therefore  be  expected  that  that 
clusters  of  grains  oriented  unfavorably  for  octahedral  slip  near  inclusions  or  voids  are  more  susceptible  to 
fatigue  crack  formation  than  inclusions  or  voids  that  are  not  located  in  such  clusters. 

A  similar  strategy  was  applied  to  investigate  the  extreme  value  distributions  of  the  driving  forces  for 
fatigue  crack  formation  in  four  different  microstructure  variants  of  duplex  Ti-6A1-4V  [142]  on 
simulations  with  fully  periodic  boundary  conditions  representative  of  subsurface  or  bulk  material 
conditions.  Additionally,  the  influence  of  a  traction  free  boundary  condition  was  also  considered  relative 
to  that  of  fully  periodic  subsurface  boundary  conditions  [147].  Multiple  SVEs  were  simulated  and  loaded 
under  strain  control  with  an  applied  strain  ratio  of  zero  (i.e.,  R£  =  0).  The  crystal  plasticity  model 
accounted  for  slip  both  the  primary  a  grains  and  the  lamellar  a+P  colonies.  In  the  primary  a  grains  that 
have  hep  lattice  structure,  four  different  families  of  slip  systems  were  explicitly  considered:  three 
(1120)(0001)  basal,  three  (1120){1010)  prismatic,  six  (1120){1011}  first  order  pyramidal  <a>,  and 
twelve  (1123){1011}  first  order  pyramidal  <a+c>  slip  systems.  Here,  slip  is  modeled  to  be  dominant  in 
the  basal  and  prismatic  slip  systems  due  to  a  relatively  low  critical  resolved  shear  stress  compared  to  the 
other  slip  systems  [81].  The  lamellar  a+P  colonies  are  homogenized  in  the  model  since  the  individual 
laths  range  in  thickness  from  hundreds  of  nanometers  to  several  microns.  A  crystallographic  Burger’s 
orientation  relation  (BOR)  is  maintained  between  the  secondary  a  and  P  laths  defined  such  that 
(0001)^  //  {101}^  and  (1120)^  //  (111)^.  The  24  possible  slip  systems  in  the  lamellar  region  include 
three  (1120)(0001)  basal,  three  (1120){1010)  prismatic,  six  (1120){1011)  first-order  pyramidal  and 
twelve  (111){110)  bcc  slip  systems.  The  bcc  slip  systems  are  transformed  into  the  hexagonal  coordinate 
system  via  the  BOR.  Hard  slip  systems  in  these  a+p  colonies  are  those  that  intersect  the  a-P  interface. 
Soft  deformation  modes  are  those  on  which  dislocations  glide  parallel  to  the  a-P  interface  or  have  parallel 
slip  planes  in  both  the  secondary  a  and  P  phases.  In  particular,  the  model  used  here  [18]  was  formulated 
to  favor  single  slip,  which  has  been  experimentally  observed  to  dominate  at  low  cyclic  strain  amplitudes. 
In  addition,  higher  critical  resolved  shear  stresses  are  required  for  slip  in  the  a+P  colonies  to  simulate 
these  slip  systems  as  is  observed  experimentally.  The  four  microstructure  variants  considered  in  this  study 
with  varying  distributions  for  the  primary  a  grain  and  a+P  colonies  are  defined  in  Table  1. 

As  shown  in  Fig.  10,  the  driving  forces  for  fatigue  crack  formation  for  the  microstructures  with  lower 
volume  fractions  of  the  primary  a  phase  and  smaller  primary  a  grain  sizes  tended  to  be  the  lowest  in 
magnitude  for  the  range  of  microstructures  considered  for  the  simulations  with  the  fully  periodic 
boundary  conditions.  Microstructure  A  was  also  simulated  with  periodic  boundary  conditions  in  two 
orthogonal  directions  with  traction  free  boundary  conditions  applied  to  two  opposing  faces  of  the  cubed 
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shaped  simulation  volume  in  the  third  orthogonal  direction.  In  this  manner,  the  simulated  volumes  were 
representative  of  a  thin  sheet  of  material  with  a  high  surface  to  volume  ratio.  It  can  be  observed  in  Fig.  1 1 
that  while  the  distributions  in  the  driving  forces  for  subsurface  fatigue  crack  formation  are  best 
characterized  by  the  extreme  value  Gumbel  distribution,  the  driving  forces  for  surface  fatigue  crack 
formation  are  best  characterized  by  the  extreme  value  Frechet  distribution.  Mathematically,  the  extreme 
value  Gumbel  distribution  describes  exponential  convergence  of  the  tails  of  the  parent  (i.e.,  initial) 
distribution;  the  extreme  value  Frechet  distribution  describes  polynomial  tails  in  the  direction  of  the 
extreme  [148].  As  the  exponential  converges  is  faster  than  a  polynomial,  the  extreme  maxima  of  the 
Frechet  distribution  are  of  greater  magnitude  than  the  extreme  maxima  of  the  Gumbel  distribution.  It 
logically  follows  that  because  the  tails  of  the  extreme  value  distributions  of  the  driving  forces  for  fatigue 
crack  formation  for  surface-dominated  cracking  are  best  described  by  polynomial  convergence  in  contrast 
to  the  extreme  value  driving  forces  for  subsurface  crack  formation  which  are  best  described  by 
exponential  decay,  the  driving  forces  for  fatigue  crack  formation  have  the  potential  to  be  highest  at  the 
surface.  However,  there  is  a  large  range  of  the  extreme  value  distributions  of  the  FIPs  for  both  surface- 
and  subsurface-dominated  fatigue  crack  formation  that  overlap.  This  overlap  in  the  extreme  value 
distribution  in  the  driving  forces  for  crack  formation  perhaps  explains  why  surface  and  subsurface  crack 
formation  are  often  both  observed  from  experiments  in  the  HCF  and  VHCF  regime  (cf.  [145]). 
Additionally,  when  the  driving  forces  for  fatigue  crack  formation  are  highest  at  the  surface  (i.e.  when  the 
extreme  value  response  at  the  surface  is  associated  with  the  extreme  maxima  of  the  extreme  value 
distribution  of  the  FIPs),  the  overall  fatigue  lives  are  predicted  to  be  shorter  than  when  fatigue  crack  form 
subsurface  where  the  extreme  maxima  of  the  extreme  value  distributions  of  the  FIPs  are  of  a  lower 
magnitude. 

The  extreme  value  marked  radial  correlation  functions  for  the  apparent  Schmid  factor  for 
Microstructure  A  of  Ti-6A1-4V  are  given  in  Figs.  11-13  for  the  both  correlations  between  the  primary  a 
grains  and  the  a+P  colonies  with  high  Schmid  factors  for  the  various  slip  systems  considered.  Here  it  can 
be  observed  that  primary  a  grains  oriented  favorably  for  basal  slip  have  a  much  higher  probability  of 
being  associated  with  the  regions  of  extreme  value  driving  forces  for  fatigue  crack  formation  than  in  the 
overall  microstructure.  Additionally,  hard-soft  grain  interactions  between  primary  a  grains  oriented 
favorably  for  basal  slip  and  primary  a  grains  oriented  favorably  for  pyramidal  <a+c>  slip  or  a+P  colonies 
oriented  favorably  for  BOR-modified  bee  slip  (i.e.,  bee  slip  systems  aligned  with  hep  systems  in 
accordance  with  the  Burger’s  orientation  relationship  (BOR))  correlated  with  increased  driving  forces  for 
fatigue  crack  formation  over  the  range  of  a+P  Ti  alloy  systems  considered.  Comparing  the  EVMCF 
generated  from  the  simulations  with  fully  periodic  boundary  conditions  versus  those  with  the  two  surfaces 
with  traction  free  boundary  conditions  demonstrated  that  the  same  coupled  attributes  likely  are  associated 
with  both  surface-  and  subsurface-dominated  failure  [147]  even  though  the  driving  forces  for  the  two 
cases  are  distributed  differently. 

The  EVMCF  provided  a  quantitative  measure  of  the  likelihood  that  correlated  microstructure  attributes 
can  be  found  in  the  neighborhood  of  the  extreme  value  response.  The  development  of  these  types  of 
statistical  tools  and  modeling  strategies  are  essential  to  describing  stochastic  processes  in  random 
heterogeneous  materials  such  as  fatigue  crack  formation  in  polycrystalline  metals.  Damage  processes  are 
particularly  sensitive  to  the  tails  of  distributions  of  the  microstructure.  By  quantifying  the  correlated 
microstructure  attributes  that  are  the  most  closely  associated  with  fatigue  damage  formation,  better 
materials  can  be  potentially  be  engineered  by  controlling  these  attributes,  as  shown  schematically  in  Fig. 
7. 


5.  Crack  Formation  and  Early  Growth  in  Notched  Specimens 

The  previous  section  focused  on  extreme  value  local  responses  induced  by  microstructure 
heterogeneity  in  the  absence  of  stress  and  strain  field  gradients.  Notches  in  aircraft  gas  turbine  engine 
components  typically  are  of  greatest  concern  with  regard  to  limiting  fatigue  life,  serving  as  preferential 
zones  for  fatigue  crack  formation  and  early  growth.  In  addition  to  the  notch,  the  effects  of  microstructure 
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heterogeneity  on  fatigue  response  remain  and  effectively  are  compounded  by  the  notch  root  gradient  field, 
giving  rise  to  notch  size  effects  that  interplay  with  both  fatigue  crack  formation  at  the  grain  scale  and 
early  growth  out  of  the  influence  of  the  notch  root  field.  Improving  the  prediction  of  variability  of  HCF 
response  of  notched  components  requires  (1)  an  understanding  of  cause  and  effect  fatigue  relationships  of 
key  microstructure  attributes,  (2)  quantification  of  notch  size  effects,  (3)  a  means  to  guide  process  route 
modification  to  enhance  fatigue  resistance,  and  (4)  experimental  validation  of  predicted  variability  of 
fatigue  life  response.  The  work  summarized  here  aims  to  address  the  first  two  items,  based  on  recent 
work  by  Musinski  and  McDowell  [149]. 

We  consider  a  supersol vus  coarse  grain  (CG)  IN  100  Ni-base  superalloy  with  an  average  grain  size  of 
34  |Lim  and  a  bimodal  distribution  of  secondary  and  tertiary  y'  precipitates,  formed  during  the  cooling 
and  aging  processing  steps  after  the  supersolvus  heat  treatment.  To  model  the  complex  mechanical 
behavior  of  the  CG  IN  100  microstructure,  a  three-dimensional,  rate-dependent,  microstructure-sensitive 
polycrystal  plasticity  model  described  in  Ref.  [141]  is  used  that  captures  the  first  order  effects  of  grain 
size  and  the  size  and  volume  fraction  of  secondary  and  tertiary  y'  precipitates.  This  polycrystal  plasticity 
model  was  calibrated  to  complex  cyclic  stress-strain  data  of  multiple  microstructure  variations  at  an 
operating  temperature  of  650°C  with  and  without  hold  times  [150]  and  was  implemented  as  a  user  defined 
material  (UMAT)  subroutine  in  ABAQUS  [107]. 

The  notched  components  in  this  study  were  double  edge-notched  specimens  with  notch  root  radius  p 
subjected  to  fully  reversed  tension-compression  fatigue  loading,  with  details  provided  in  Ref.  [149].  The 
model  employed  a  finer  mesh  size  of  4  |Lim  in  a  localized  region  near  the  notch  root  (0.15/7  away  from 

notch  root)  and  contained  a  graded  mesh  that  gradually  increased  to  the  CG  size  (34  |Lim)  at  the  outside 
barrier  of  the  poly  crystal  plasticity  zone.  The  poly  crystalline  grain  structure  within  the  notch  root  region 
was  constructed  using  a  random  sequential  adsorption  algorithm  similar  to  that  previously  described  in 
Refs.  [123-124,129].  The  values  of  p  -  -0.1  and  a  -  0.4  were  chosen  for  the  target  lognormal  grain  size 

(mean  grain  size  =  34  |um)  distribution  function,  f(x\  //,  <J)  — - 

X<TV  2  71 

previous  publications  of  fine  grain  IN100  grain  size  distributions  [83,123-124,  141,151].  An  example  of 
the  polycrystalline  grain  structure  is  shown  in  Fig.  14(a),  where  each  grain  is  represented  by  a  different 
color  to  visualize  the  grain  structure.  The  cumulative  Von  Mises  effective  plastic  strain  within  the  notch 
root  region  at  maximum  tensile  load  during  the  third  fatigue  cycle  is  shown  in  Fig.  14(b).  In  this  Figure, 

the  values  of  effective  plastic  strain  below  a  value  oi£p  =  1CT7  were  dropped  from  the  contour  plots  to 
display  the  heterogeneity  of  plastic  slip  in  the  microstructure.  As  expected,  the  effective  plastic  strain  is 
most  significant  at  certain  grain  boundaries,  which  results  from  the  incompatibility  of  deformation 
between  grains  of  different  size  and  orientation.  Plotted  in  Fig.  14(c)  is  the  distribution  of  maximum  FIP 
with  respect  to  x-distance  from  the  notch  root  for  five  different  notch  root  sizes  at  strain  amplitude 
£a  =  £y  ,  where  £  is  the  yield  strain. 

Since  the  maximum  FIPfs  =  AFmax(x)  (Eq.  (1))  as  a  function  of  x-distance  from  the  notch  root  is 


(ln(x)-//)2 

2c-2 


,  based  on 


relatively  linear  on  a  log-linear  plot  (Fig.  14(c)),  AFmax(x)  was  parameterized  as  an  exponential  decay 


function,  Arm^(x)  =  ArmOT(0)exp 


V  Ld  J 


where  <^is  a  constant  and  L  ,  is  the  transition  crack  length. 


This  log-linear  best  fit  was  used  to  estimate  MSC  growth  within  the  notch  root  field,  via 
=  Afs  Ty  Ar max(x)  a  ,  where  T^is  the  critical  resolved  shear  stress  and  a  is  the  crack  length.  This 


relation  was  used  to  construct  a  crack  length  (an)  versus  number  of  cycles  ( N )  relationship,  and  to 
assess  the  probability  of  failure  at  any  given  number  of  cycles  via  the  cumulative  distribution  function 
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CDF(Sa-,Ld,N)  =  1-exp 


;  here, 


77  =  —  ln(0.5)  was  used  to  normalize  the  CDF  so  that 


CDF=0.5  (50%  probability)  when  the  crack  extends  to  the  transition  crack  length,  Ld  .  This  method  of 

determining  the  CDF  function  is  summarized  in  Fig.  15. 

Figure  16  compares  strain-life  data  and  simulation  results  for  smooth  and  notched  specimens.  Details 
of  the  smooth  specimen  computational  fits  to  experimental  data  [152-154]  in  Ref.  [149].  As  seen  in  the 
Fig.  16,  a  notch  effectively  knocks  the  strain-life  curve  downward  in  the  HCF  and  VHCF  regimes.  It  is 
interesting  to  note  in  Fig.  16  that  the  poly  crystal  plasticity  model  predicts  a  shallow  slope  in  the  VHCF 
regime  for  smooth  specimens  based  on  calibration  to  experimental  data  in  LCF  and  LCF-HCF  transition 
fatigue  regimes.  This  is  a  consequence  of  constrained  microplasticity  at  low  strain  amplitudes  [149]. 

The  results  for  the  CDF  are  shown  in  Fig.  17  for  a  strain  amplitude  of  £a  =  0.6^ .  Each  blue  dot  in 

Fig.  17  represents  the  CDF  evaluated  at  a  given  number  of  cycles,  N  ,  for  an  individual  instantiation  of  a 
notched  specimen  with  a  random  microstructure.  Therefore,  since  20  different  instantiations  were 
simulated  per  notch  root  radius  at  this  strain  amplitude,  there  are  20  blue  dots  for  a  given  x-axis  value  of 
life,  N  .  The  transition  from  MSC  growth  to  LEFM-correlated  long  crack  growth  is  denoted  by  Pf 

values  above  0.5.  The  black  symbols  in  Fig.  17  represent  the  mean  Pf  value  for  all  20  instantiations  at  a 

given  number  of  cycles.  As  can  be  seen  from  these  Figures,  the  benefit  of  using  this  approach  is  that  the 
cumulative  distribution  function  can  be  calculated  for  any  number  of  cycles  and  any  probability  of  failure. 
This  approach  can  be  applied  to  LCF,  HCF  and  transition  fatigue  regimes  to  compute  the  total  component 
failure.  Once  the  crack  forms  to  the  transition  length,  LEFM  can  be  used  to  assess  crack  growth  to 
specimen  failure.  This  framework  predicts  that  larger  notch  sizes  will  tend  to  fail  before  smaller  notch 
sizes,  consistent  with  general  experimental  trends,  in  addition  to  microstructure  size  effects  on  crack 
formation  and  early  growth  within  the  influence  of  the  notch  root.  These  computationally  estimated 
cumulative  distribution  functions  need  to  be  further  compared  with  (limited)  experiments  for  validation 
and  further  refinement. 


6.  Closure 

Recent  progress  in  mesoscopic  modeling  and  simulation  that  employs  polycrystal  plasticity  as  the 
primary  computational  tool  has  been  examined  for  applications  in  fatigue  of  Ti  and  Ni-base  superalloys, 
with  emphasis  on  the  HCF  and  VHCF  regimes  which  are  dominated  by  crack  formation  in  individual 
grains/phases  and  early  growth  in  the  neighboring  microstructure.  Following  a  brief  review  of  the 
history  of  mesoscopic  simulations  that  seek  to  quantify  the  role  of  material  heterogeneity  on  fatigue 
response,  a  simple  methodology  is  outlined  based  on  fatigue  indicator  parameters.  It  is  argued  that  the 
cyclic  crack  tip  displacement  serves  as  a  useful  construct  as  a  driving  force  for  crystallographic  fatigue 
crack  growth  that  bridges  scales  from  individual  grains  to  crack  growth  over  larger  numbers  of  grains 
based  on  EPFM  concepts.  It  is  shown  by  direct  simulation  of  cyclic  loading  of  stationary  cracks  in  single 
crystals  that  the  Fatemi-Socie  FIP  directly  relates  to  the  ACTD.  Variability  in  HCF  and  VHCF  is 
computationally  projected  by  using  multiple  statistical  volume  elements  and  considering  the  extreme 
value  distributions  built  up  from  these  simulations.  The  concepts  of  marked  correlation  functions  and  are 
reviewed  as  a  means  to  quantify  the  role  of  multiple  microstructure  attributes  that  couple  to  produce  the 
extreme  value  FIPs  in  the  HCF  regime.  An  algorithm  is  outlined  for  estimating  the  cumulative 
probability  distribution  of  cycles  for  crack  formation  and  growth  from  notches  in  HCF  and  VHCF. 

This  class  of  approaches  is  promising,  but  much  work  remains  to  be  done  in  terms  of  validation  for 
more  extensive  datasets  and  materials,  as  well  as  extensions  that  consider  processes  of  fatigue  crack 
nucleation  in  more  detailed  and  complete  fashion,  including  the  role  of  grain/phase  boundaries,  and 
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methods  to  faithfully  simulate  growth  of  microstructurally  small  cracks  in  fatigue  through  3D  ensembles 
of  grains.  Relatively  little  attention  has  been  devoted  to  these  issues. 
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Tables 


Table  1.  Microstructures  for  Ti-6A1-4V  along  with  the  assumed  log-normal  fits  for  the  grain  volume 
distributions. 


Duplex  Ti-6A1-4V  Microstructures 

Assumed  Mean  and  St.  Dev.  for 
Grain  Size  Distributions 
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P  Size  (pm) 

Primary  a 
Size  (pm) 

Vol  % 
Primary  a 

(pm) 

<5 

(pm) 

ii 

(pm) 

a 

(pm) 

A 

Fine  bi-modal  low  a 

50 

10-50 

30% 

50 

5 

25 

10 

B 

Fine  bi-modal  high  a 

50 

70% 

50 

5 

25 

10 

C 

Coarse  bi-modal  low  a 
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Figures 


Crack  formation  due  to  intense 
shear  along  the  slip  band  of  Ti-6AI- 
4 V  [75]. 


Slip  band  impingement  on  grain 
boundary  of  polycrystalline  nickel 
[76]. 


Fig.  1.  Candidate  fatigue  indicator  parameters  (FIPs)  with  combined  plastic  shear  strain  and  normal  stress 
effects  on  candidate  planes  for  crack  nucleation  and  early  growth  along  slip  bands  or  grain  boundaries  for 
Ti-6A1  [75]  and  Ni  [76]. 
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Simulations 


Distribution  of  FIPs 


Digital  SVEs 

t 

Robust  Materials 
(Morphology)  Design 


Experimental 

Calibration/Validation 


Rank-Order  Range  of 
Microstructures  (for 
max./min.  life  or  min. 
variability) 


Assess  Fatigue 
Variability  via  FIPs 


Fig.  2.  Overall  concept  of  using  simulations  from  multiple  SVEs  instantiated  using  spatial  statistics  of 
representative  microstructures  to  compute  the  distribution  of  FIPs  and  estimate  the  variability  of  fatigue 
response,  using  this  information  to  rank  order  microstructures  and  to  provide  feedback  to  process  route 
and  structure  to  improve  material  performance  [99]. 


Fig.  3.  Constant  life  contour  plots  in  Tp  -plane  (based  on  plastic  strain)  for  two  different  LCF  lives  based 
on  the  Fatemi-Socie  parameter  (reproduced  from  Ref.  [63]). 
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Fig.  4.  Cracked  single  crystal  employed  to  evaluate  the  fine  scale  FIP  and  ACTD  measured  between 
nodes  in  red  at  the  mid-thickness  for  cyclically  loaded  stationary  cracks.  The  crystallographic  planes  are 
oriented  so  that  the  crack  plane  is  an  octahedral  slip  plane  with  a  slip  direction  along  the  X-axis. 


Fig.  5.  Shear  (left)  and  mixed  mode  (right)  loading  cases  employing  the  same  magnitude  of  the 
displacement  vector  (D). 
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Fig.  6.  Comparison  of  computed  fine  scale  FIP  versus  ACTD  for  (left)  shear  and  (right)  mixed  mode 
loading  for  Cu  (blue  symbols)  and  RR1000  Ni-base  superalloy  (black  symbols)  with  small  cracks  under 
cyclic  loading;  solid  lines  have  unity  slope,  indicating  a  linear  relation. 


1.  Generate  multiple  SVEs 
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distributions  of  key 
microstructure  attributes 
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2.  Identify  EV  response  of  SVEs 
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candidates  for 
experimental 
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5.  Identify  extreme  value  correlated 
attributes  key  to  response  and  rank 
microstructures 


3.  Characterize  EV 
distributions  of  key  response 
parameters 

4.  Characterize  correlated 
microstructure  attributes 
coincident  with  the  EV 
response  (EV  marked 
correlation  functions) 


Fig.  7.  Schematic  of  a  microstructure-sensitive  probabilistic  framework  of  extreme  value  type  in  which 
simulations  of  many  SVE  instantiations  provides  feedback  to  processing  and  materials  develop 
microstructures  with  improved  performance/variability  in  HCF  [99]. 
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Extreme  Value  Fatemi-Socie  FIP 


Fig.  8.  For  IN100,  the  grain  averaged  extreme  value  Fatemi-Socie  (FS)  FIP  for  the  SVEs  cycled  at  0.3%, 
0.5%,  and  0.7%  maximum  strain  are  plotted  on  Gumbel  probability  scales.  Note  that  the  equations  given 
for  the  least  squares  linear  regression  are  such  that  y  =  ln(l/ln(l/p))  and  x  =  extreme  value  FIP.  The 
extreme  value  FIP  for  each  SVE  was  selected  to  as  the  maximum  of  those  calculated  over  each  grain. 
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(a) 


r/d 


(b) 

Fig.  9.  For  IN100  the  (a)  complete  radial  and  extreme  value  marked  radial  distribution  for  apparent 
Schmid  factors  m8  =0.45  to  0.5  (cube  slip)  and  m8a  =0.45  to  0.5  (cube  slip)  and  (b)  complete  radial  and 

extreme  value  marked  radial  distribution  of  apparent  Schmid  factors  m8  =0.45  to  0.5  (cube  slip)  and  m8 
=0.45  to  0.5  (octahedral  slip)  are  plotted  for  the  100  SVEs  subjected  to  0.5%  strain.  The  distance  r  that 
separates  the  two  orientations  is  normalized  against  the  average  cube  root  grain  size  d  (i.e.,  20  |Ltm). 
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Fatemi-Socie  FIP 


Fig.  10.  The  extreme  value  Fatemi-Socie  FIPfs  calculated  over  a  cube  shaped  averaging  volume  with 
equivalent  grain  size  of  32  |Lim  for  the  four  analyzed  Ti-6A1-4V  microstructures  plotted  on  Gumbel  (Type 
I)  probability  scales. 
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(a) 


(b) 


Fig.  11.  The  extreme  value  Fatemi-Socie  FIPfs  for  Ti-6A1-4V  calculated  over  a  cube  shaped  averaging 
volume  with  equivalent  grain  size  of  0.032  |Lim  for  the  instantiated  SVEs  (R8=0  and  maximum  applied 
strain  of  0.006)  with  periodic  BCs  and  with  free  BCs  plotted  on  (a)  Gumbel  (Type  I),  (b)  Frechet  (Type 
II)  probability  scales.  The  equations  given  for  the  least  squares  linear  regression  are  such  that  y  = 
ln(l/ln(l/p))  and  x  =  FIP  for  the  Gumbel  probability  plot  and  x  =  In(FIP)  for  the  Frechet  probability 
plot,  respectively. 
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Fig.  12.  Radial  correlation  functions  for  the  complete  Ti-6A1-4V  microstructure  compared  with  EVMCF 
describing  the  correlation  between  the  apparent  Schmid  factor  m8  for  basal  slip  between  0.45  and  0.5  for 

the  primary  a  phase  and  the  Schmid  factor  m8a  for:  (a)  basal  slip  between  0.45  and  0.5,  (b)  prismatic  slip 

between  0.45  and  0.5,  (c)  pyramidal  (a)  slip  between  0.45  and  0.5,  and  (d)  pyramidal  (a  +  c)  slip 

between  0.45  and  0.5  for  the  primary  a  phase  for  the  instantiated  SVEs  with  fully  periodic  boundary 
conditions. 
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Fig.  13.  Radial  correlation  functions  for  the  complete  microstracture  compared  with  EVMCF  describing 
the  correlation  between  the  apparent  Schmid  factor  mga  for  basal  slip  between  0.45  and  0.5  for  the 

primary  a  phase  and  the  apparent  Schmid  factor  m8a  for:  (a)  basal  slip  between  0.45  and  0.5,  (b) 

prismatic  slip  between  0.45  and  0.5,  (c)  pyramidal  (a}  slip  between  0.45  and  0.5,  and  (d)  (l  1  l){l  10} 

bcc  slip  (transformed  into  the  hexagonal  coordinate  system  via  the  BOR)  between  0.45  and  0.5  for  the 
a+P  colony  phase  for  the  instantiated  SVEs  with  fully  periodic  boundary  conditions. 
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(a)  Polycrystalline  grain  structure. 


(b)  Effective  plastic  strain. 
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(c)  FIP  distribution  as  a  function  of  X-distance  from  notch  root. 


Fig.  14.  (a)  Polycrystalline  grain  structure,  (b)  Von  Mises  effective  plastic  strain  for  a  coarse  grain  IN100 
notched  specimen  at  maximum  tensile  load  during  the  third  fatigue  cycle,  p  =  0.6  mm,  £a  -  0.6£y ,  and 
(c)  maximum  FIP  distribution  versus  x-distance  from  notch  root  for  five  different  notch  root  sizes  for  the 
coarse  grain  IN  100  microstructure  at  a  strain  amplitude  £a  =  £y  . 
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Fig.  15.  Overview  of  method  for  assessing  CDF  for  the  transition  crack  length  approach  based  on 
uncracked  FIP  distribution  (details  in  Ref.  [149]). 
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Figure  16.  Comparison  of  coase  grain  IN100  smooth  and  notched  computational  strain  life  data  versus 
HIP  Astroloy  (Cowles  et  al.,  1978,  1980)  and  N18  (Bathias  and  Paris,  2010)  smooth  specimen  strain-life 
experimental  data. 
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Fig.  17.  Coarse  grain  IN100  cumulative  distribution  function  for  various  notch  root  radii  ( £a  =  0.6  e  , 
oa  =  450 ,  Re  =  -1 ,  T  =  650°C ). 
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